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The fomalism is developed to express nucleon matrix elements of the electromagnetic current in terms of form 
factors consistent with the translational, rotational, and parity symmetries of a cubic lattice. We calculate the 
number of these form factors and show how appropriate linear combinations approach the continuum limit. 


Electroweak form factors provide a precise ex¬ 
perimental probe of the quark substructure of 
hadrons, so it is of great interest to calculate 
them in QCD from first principles. Since the 
only known way to solve QCD is numerical so¬ 
lution on a discrete lattice, this work addresses 
the decomposition of matrix elements of the elec¬ 
tromagnetic current into form factors consistent 
with cubic lattice symmetry and the relation of 
these form form factors to the familiar ones asso¬ 
ciated with continuum Lorentz symmetry. 

The breaking of continuum symmetries leads, 
in general, to more form factors on the lattice 
than in the continuum. In a practical Monte 
Carlo calculation on the lattice, it is advanta¬ 
geous to calculate these lattice form factors since 
all momenta related by a lattice symmetry can 
be combined to obtain the highest statistical ac¬ 
curacy. Therefore, it is necessary to determine 
the appropriate combination of these form fac¬ 
tors that approaches the physical form factors in 
the continuum limit. 

In this work, we consider the electromagnetic 
form factors of the nucleon. The matrix element 
of the electromagnetic current in the continuum 
is 

{p',s'\J^\p,s) = 

Hp',s') u{p,s) 
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where p and s {p' and s') are the initial (final) 
momentum and spin z-component and is the 
electromagnetic current. This is the most gen¬ 
eral form consistent with continuum space-time 
symmetries and conservation of the electromag¬ 
netic current. We seek an analogous result on the 
lattice. 

It is easiest to consider the following lattice ma¬ 
trix element between spin and helicity states: 

{6,s\Jf,\p,h). 

The final state has zero momentum and spin z- 
component s. The initial state has an arbitrary 
non-zero momentum p and helicity h. The fact 
that helicity states do not mix under rotations 
simplifies our analysis. We now need to construct 
the most general form of this matrix element con¬ 
sistent with lattice spatial symmetries. 

To construct lattice form factors, we decom¬ 
pose the non-zero momentum state into irre¬ 
ducible representations of the lattice rotation 
group and apply the Wigner-Eckart theorem to 
write the matrix element in terms of a set of re¬ 
duced matrix elements. These reduced matrix 
elements are the lattice form factors. Since the 
final state is already irreducible and the current 
simply consists of two irreducible representations, 
Jo and Ji, 2 , 3 , we only need the decomposition of 
the non-zero momentum states and the Clebsch- 
Gordan coefficients. 

To start, we consider the orbit of momentum p 
under the lattice rotation group, that is, all lattice 
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Table 1 

For each type of momentum class, where N > 
M > L > 0, we list the order of the class, the 
number of ^ and | representations in the decom¬ 
position, the number of G and F lattice form fac¬ 
tors, and the total number of lattice form factors. 


class 

order 



G 

F 

G + F 

[(0,0,0)] 

1 

1 

0 

1 

1 

2 

[(iV,0,0)] 

6 

1 

1 

1 

2 

3 

[iN,N,N)] 

8 

1 

1 

1 

2 

3 

[{N,N,0)] 

12 

1 

2 

1 

3 

4 

[(7V,M,0)] 

24 

2 

4 

2 

6 

8 

[{N,M,M)] 

24 

2 

4 

2 

6 

8 

m,N,M)] 

24 

2 

4 

2 

6 

8 

[{N,M,L)] 

24 

2 

4 

2 

6 

8 

[iN,M,-L)] 

24 

2 

4 

2 

6 

8 


rotations of p. This set of states forms a basis of a 
generally reducible representation in the Hilbert 
space. We block diagonalize this representation, 
producing a change of basis from linear to angular 
momentum states that may be written as follows: 

= 

8 riR dg 

EEE Ac,h{p\ R, a, s)\R, a, 

R—1a—1 s—1 

where R labels the eight irreducible, double¬ 
valued representations of the universal cover of 
the lattice rotation group (Ref. 0.i ), nn is the 
number of occurrences of R, a labels each occur¬ 
rence, du is the dimension of R, s labels the basis 
states, C denotes a momentum class from which 
p takes its values, and h is the helicity. All the 
momentum classes (except for irrelevant special 
cases at the edge of the Brillouin zone) can be 
catalogued into the eight cases shown in Table 1, 
of which only hve have essentially different A’s. 
Note that all momenta are lattice momenta, and 
physical momenta require a factor of tt/L. 

The general results for matrix elements may be 
written as follows: 


{0,s\Jo\p,h)(. = 

ni 

OC 

{6,s\Ji\p,h)(. = 

n 1 

Oi 

ns 

E ^c,h{p'^ 3/2, a, s, 

Oi 

Bc,hip;R,a,sG) = 

dR 

E ^c,h{p\ R, a, s')C'(l/2, s, 1, i] R, s) 

s' 

where C(l/2, s, 1, z; i?, s') denotes the Clebsch- 
Gordan coefficients coupling ^ 0 1 to i?. The 
states (0, s|Jo transform as the i representation 
and hence have overlap only with the ^ compo¬ 
nent of \p,h)Q. Therefore for each occurrence of 
the ^ representation in \p, h)^ there is an unde¬ 
termined reduced matrix element, G. The states 
(0, s| Ji transform as the ^ 0 1 = ^ 0 | represen¬ 
tation and hence overlap only with the ^ and | 
pieces of |p, h)Q. Similarly for each occurrence of 
the ^ or I representation in the decomposition 

there is an undetermined matrix element, Fi or 

2 

F3. 
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To illustrate the above results, we discuss two 
examples. The first treats momenta for which 
the continuum form of the Jq rnatrix element is 
preserved on the lattice. This occurs for several 
types of momentum classes: [(0,0,0)], [(A^, 0,0)], 
[(A^, A^, A^)], and [(A^, A^, 0)]. The result for mo¬ 
menta of the type [(A^, 0,0)] is 

{0,s\Jo\p,h) =-^Xsihp)G[^^ 

where Xsip) denotes a standard spin one-half 
spinor polarized in the p direction and has 
the following continuum limit 


/^lat ^ ^ fo j /^cont 

Note that there are arbitrarily large lattice mo¬ 
menta for which the continuum form is retained. 
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The second example treats momenta for which 
the Jo matrix element has an additional form fac¬ 
tor. This occurs for the remaining types of mo¬ 
mentum classes, and we consider momenta of the 
type [{N, M, 0)]. Since there are two independent 
terms G^, we have taken a linear combination 
such that one term is proportional to the spinor 
Xs(hp), with the result 


(0,s|Jo|p, h) = 

^Xs{hp)G^r + ^Xs{-hp)G'r- 


Hence, the two form factors have the continuum 
limit 


^lat 




G 


cont 

e 


and we have succeeded in identifying the lin¬ 
ear combination that corresponds to continuum 
physics. 

Note that on the cubic lattice, in addition to 
the ^ representation corresponding to the physi¬ 
cal angular momentum i state in the continuum, 
there exists a second i representation originating 
from higher angular momentum continuum states 
that are imperfectly represented on the cubic lat¬ 
tice. Since there is no physical matrix element 
in the continuum connecting the ^ state to these 
higher angular momentum states, it follows that 
the coefficient 6 * 2 ^* must vanish in the continuum 
limit. 

The relation between behavior on the lattice 
and in the continuum can also be clarified by con¬ 
sidering the little group for each non-zero momen¬ 
tum. The little group, or stabilizer, is the sub¬ 
group of lattice rotations which leave the momen¬ 
tum invariant. Each irreducible representation of 
the little group leads to an allowed helicity, as 
summarized in Table 2. Just as the number of dis¬ 
tinct spins is restricted on the lattice, the number 
of allowed helicities is also limited and dependent 
upon momentum. In the continuum G^ 2 ^ would 
mediate an ^ ^ ^ transition, but on the lattice, 

for this type of momentum class, — ^ is isomor¬ 
phic to i, thus the lattice does not have enough 
constraints to exclude G^f^. 


Table 2 

For each type of non-zero momentum class we 
list: the little group, the allowed helicities, and 
the helicity range. 


class LG helicities Ah 

[(iV,0,0)] Zg -1,-i, 0, i, 1, |, 2 r 

[{N,N,N)] Ze -1,-i, 0, i, 1, I 3 

[{N,N,0)] Z4 -iO, i, 1 2 

[(iV,M,0)] Z 2 0, i 1 

[{N, M, M)] Z 2 0, I 1 

[(A,iV,M)] Z 2 0 , I 1 

[{N,M,L)] Z 2 0, I 1 

[{N,M,-L)] Z2 0,1 _^ 


The previous example demonstrates a general 
feature: the order of the little group determines 
the extent to which the continuum form is iden¬ 
tically adhered to by the lattice. A larger little 
group ensures sufficient symmetry to enforce the 
continuum results, whereas a smaller little group 
lacks the constraints to prevent additional form 
factors. 

To summarize, the lattice breaks continuum 
symmetries, and hence the number of lattice form 
factors increases. The remaining symmetries vary 
with momentum, consequently, there are mo¬ 
menta which admit one to six extra form factors. 
We can understand the presence of these addi¬ 
tional form factors in terms of the degraded range 
of helicities for each momentum. The explicit 
decomposition of lattice matrix elements allows 
us to distinguish these unphysical lattice form 
factors. Our future work will apply this analy¬ 
sis to the SESAM configurations and compare it 
with the conventional analysis using continuum 
expressions. 
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